A = rand(4);

exact = expm(A);

% Error in taylorsimple

T = taylorsimple(A);
taylorfout = norm(exact-T)

% Error in eigsimple

S = eigsimple(A);
simplefout = norm(exact-S)

% Error in eigqr

Q = eigqr(A);
qrfout = norm(exact-Q)